Jaynes-Cummings dynamics with a matter wave oscillator 
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We propose to subject two Bose-Einstein condensates to a periodic potential, so that one con- 
densate undergoes the Mott insulator transition to a state with precisely one atom per lattice site. 
We show that photoassociation of heteronuclear molecules within each lattice site is described by 
the quantum optical Jaynes-Cummings Hamiltonian. In analogy with studies of this Hamiltonian 
with cavity fields and trapped ions, we are thus able to engineer quantum optical states of atomic 
matter wave fields and we are able to reconstruct these states by quantum state tomography. 
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The Jaynes-Cummings model was introduced to de- 
scribe the resonant interaction between an atom and a 
single mode of the quantized radiation field fl , and with 
incorporation of pumping and relaxation mechanisms, 
it constitutes the corner stone for the quantum treat- 
ment of maser and laser action |2| . The unitary Jaynes- 
Cummings dynamics has been succesfully demonstrated 
in the micro maser in which atoms in Rydberg ex- 
cited states are injected into a high-Q cavity with a single 
mode which is resonant with an atomic transition. This 
system has been used to demonstrate fundamental as- 
pects of quantum theory such as the oscillatory exchange 
of energy between the atom and the field, sub-poissonian 
photon statistics, entanglement, and decoherence [Q, |). 
Another implementation of the Jaynes-Cummings Hamil- 
tonian is in the description of the quantized motion of an 
ion in a harmonic trap ||, for which a classical laser 
field drives transitions in which the ion is excited and 
simultaneously one quantum of energy is extracted from 
the motional state. The ion system is characterized by 
very long lifetimes of the oscillator state, and hence the 
Jaynes-Cummings dynamics has been studied in great 
detail, and numerous examples of state synthesis and de- 
tection have been presented @, |§| . 

In a quantum optical perspective, what characterizes 
a Bosc-Einstein condensate is its single-mode character 
and coherence properties ||. The bosonic commutator 
relation applies for atoms as well as for photons and 
phonons, and if one neglects interactions, many features 
of single mode fields are thus common for atoms, light 
and harmonic motion. This has of course led to various 
proposals for generation of atomic states with proper- 
ties similar to the ones of light, e.g., phase and number 
squeezed states and Schrodinger cat states. By clever use 
of the atomic interactions as an adjustable non-linear el- 
ement, strong quantum effects are foreseen ||l0| , and in 
combination with the long storage times for atoms, these 
effects are potentially much more useful, for atoms than 
for light, and indispensable, e.g., for the improvement of 
atomic clocks [ jTl| . 

In this Letter we present a quantum optical analysis 



of the dynamics of a pair of atomic condensates which 
have been trapped in a periodic potential, so that one 
species has experienced the Mott insulator transition to 
a state with exactly one atom per lattice site Jl^, 
The other condensate is assumed to have a different in- 
teraction strength, and hence it may still be in the su- 
pcrfluid phase with Poissonian number fluctuations on 
each site, or it may have reduced fluctuations due to the 
on site repulsion. A detailed presentation of this process 
was given recently pjfl . When the optical lattice poten- 
tial is raised the atomic wave functions become more lo- 
calized in the potential minima, and thus the collisional 
interaction between atoms increases and the tunneling 
rate between wells decreases. From the time dependen- 
dent shape of the potential for the atoms one determines 
the on-site interaction strengths U ai Ub and U a b between 
atoms of the same and different species A and B and the 
interwell tunneling couplings J a , Jb, but for simplicity 
we shall just assume the values shown in the upper panel 
of Fig. 1. The ratio of approximately 0.6 between U a and 
Ub corresponds to the ratio between the free space s-wave 
scattering lengths of Rb and K atoms. The tunneling 
coupling of both species are taken to follow the same 
lower curve. These curves are in qualitative agreement 
with the dependence of the interaction strength and tun- 
nel coupling for a single species evaluated in [jlB), where 
the energy and frequency unit Jo, on the order of a tenth 
of the atomic recoil shift, is introduced (Ti = 1). Using 
a symmetry breaking Gutzwillcr ansatz, we assume that 
the full many body system is in a product state over 
all sites, where each site is described by the state vec- 
tor = Y, nA ,n B a nA ,n B \nA,n B ), with the number 
of atoms of species i = A, B. We have solved the self- 
consistent equations for the amplitudes a„ Ai „ B , and we 
have verified that the two species undergo the Mott insu- 
lator transition at different times, see also |Q. The lower 
panel of Fig. 1. shows typical results of a numerical so- 
lution of the problem. The solid lines show the variance 
of the populations (the mean values are unchanged and 
equal to the initial variance where we assume Poisson 
statistics). The dashed curves are obtained by neglect- 
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ing the cross coupling between the two species, and as 
in juj we verify that this coupling has little effect on 
the population statistics of the two species in the lattice 
wells. 




FIG. 1: Mott insulator transition for two atomic species. The 
upper panel shows the time dependent coupling constants 
and the tunneling coupling, all in units of the constant Jo. 
The lower panel shows the variance of populations of the two 
atomic species. The dashed curves are obtained by neglect of 
the cross coupling term U a b- ■ 

Numerous studies have addressed the prospects of 
forming molecules from ultra-cold atoms both theoret- 
ically and experimentally It has been suggested to 
trap the atoms in lattices with low occupation number, 
and either make explicit use of the Mott insulator transi- 
tion 14 15f, or simply work with low mean occupation of 
each site 17 to form the molecules in an optimally con- 
trolled manner. We suggest to build on these proposals 
and to drive the photoassociation process on each lattice 
site to form a heteronuclear molecule with a single atom 
from each of the two species. We are more interested in 
the process than in the final molecular state, and hence 
a sufficiently long lived excited vibrational state will suf- 
fice as molecular state of the process. The coherent pho- 
toassociation process is described by the effective, second 
quantized, Hamitonian 
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where a,i is the annihilation operator of atoms and 
molecules (i = A,B,M). When restricted to the states 
with initially precisely one B atom, we can rewrite this 
Hamiltonian as 
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where the Pauli spin lowering operator a = a B a,M 
represents the two-state transition operator 
|1b) 0m)(0b, 1m|j (<7 turns an AB molecule into 



a B atom, i.e. , it does not conserve the number of A 
atoms, and it appears only in combination with the 
creation operator a' A ). The Hamiltonian couples pairs of 
levels 



\n A , 1 B , Om) <-> \(ii-1)a, Ob 1m), 



(3) 



and it thus follows that there will be an oscillatory ex- 
change of the single B atom between being part of the 
molecule and being an independent atom, in complete 
analogy with the Jaynes-Cummings transfer of excitation 
between a two-level system and an oscillator mode. 

This oscillation will undergo collapses and revivals, be- 
cause the coupling strength betwen the two mentioned 
states is proportional to ^/ua, and hence the different 
number state components of A atoms will give away an 
atom to the association process at different frequencies. 
An experiment that counts the total number of B atoms 
as a function of duration of the photoassication process, 
can thus teach us about the number statistics of the A 
atoms! 

So far, the Mott insulator transition has been demon- 
strated experimentally by the loss of interference between 
different lattice sites, and by a gap in the excitation spec- 
trum [ fl3|| . The Jaynes-Cummings dynamics presents a 
unique application of a quantum optical method for the 
study of this many-body problem. We suggest to turn 
on slowly the lattice potential, so that the B atoms un- 
dergo the transition to the Mott phase at a time where 
the A atoms are still in an unknown state. At this time 
we hold the lattice potential and we turn on the pho- 
toassociation process, i.e., the dynamics driven by the 
Jaynes-Cummings Hamiltonian (Q). The probability to 
find a B atom on any lattice site is given by the square 
of the Rabi oscillation amplitudes, weighted by the ua 
probability distribution p{ua)'- 



(n B ) 



p(n A ) cos 2 (xVnAt) 



(4) 



This function is illustrated in Fig. 2. which shows the 
mean value of ub as a function of time in case of cou- 
pling to a Poissonian distribution with on average two 
A atoms per lattice site (upper panel), a sub-Poissonian 
distribution of A atoms achieved at Jot = —50, cf. Fig 
1, during application of the periodic potential (middle 
panel), and an ua = 2 (Mott insulator) number state 
(lower panel). For time sequences like the ones in Fig. 2, 
one can use Eq.(^) to obtain a ht, which produces the 
complete p{ua) distribution. By counting the total num- 
ber of B atoms, without requiring the experimental abil- 
ity to detect atoms with high efficiency on a microscopic 
spatial scale, we thus obtain the precise number statis- 
tics at each site of the A atoms. A similar procedure was 
initially applied for trapped ions Q, where the internal 
state of the ion was used to monitor the population of 
the harmonic oscillator levels of motion. Many repeated 
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experiments on the same system were necessary to probe 
the state of the ion at different durations t of the dynam- 
ics and with good precision. We also have to operate the 
process many times with different durations of the pho- 
toassociation process, but since the experiment 'repeats 
itself in parallel in every lattice site, determination of 
the total number of B atoms in the entire system suffices 
to give the mean value of tlb- To take into account both 
deviations from the Lamb-Dicke limit and decoherence, 
in Ref. 0| the ideal Jaynes-Cummings dynamics ([l|) was 
replaced by a more general formula 

(n B ) = ^ p{n A ) cos 2 {n nA t)e-"< n a 1 (5) 

with fitting parameters p(tia), ^n A and "f nA - Similar 
extra handles on the theory may be useful in our matter- 
wave analog where, e.g., atom-atom interactions may 
cause an n^-dependence of the atomic wave function and 
hence of the photo-association coupling rate, and where 
various mechanisms, such as inhomogeneities over the 
size of the lattice, may cause decoherence. 
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FIG. 2: Time dependence of the mean number of B atoms per 
lattice site during photoassociation/-dissociation. The upper 
panel shows the results for a Poissonian distribution with an 
average of two A atoms per site; the middle panel shows the 
results for the distribution of A atoms at Jot = —50 in Fig. 
1, where the Mott insulator transition has not yet taken fully 
place for the A atoms; and, the lower panel shows the results 
for an ua = 2 number state. 

In case one can only effectively control the Mott- 
insulator transition of one species (the A atoms), we note 
that for a sufficiently low density of B atoms, every lattice 
site will be occupied by a single or none of these atoms, 
and hence the mean value of ns will ocillate in the same 
way as in (BUB)) but multiplied with the mean number of 
B atoms per site. 

As mentioned in the introduction of the Letter, the 
Jaynes-Cummings Hamiltonian has been applied exten- 
sively in quantum optics, and the range of interest- 
ing phenomena studied is very wide. Let us turn to 



the possibilities for synthesis and studies of particu- 
lar quantum states of matter by means of the effective 
Jaynes-Cummings Hamiltonian. Studies of the collapse 
and revival of the diffraction pattern of lattice trapped 
atoms have already been made Jll| |2(J, and we pro- 
pose to use the Jaynes-Cummings dynamics to produce a 
Schrodinger cat like state of the A atoms, and to probe its 
coherence properties. The photoassociation process can 
be driven off resonance with a detuning S in which case 
the two-state transition (^) does not take place, but the 
coupling (H) perturbs the energy of the state \tia, Is) 
by the amount AE nA = \x\ 2 ua/S. If the A atoms are 
in a coherent state |a), this state will undergo a simple 
phase rotation, | a) — > \ae 1 ^) due to this coupling. Imag- 
ine now that the B atom has a level structure, so that 
it can be initally prepared in a superpostion of the ac- 
tive state \B) which experiences the coupling discussed 
above, and a passive state \B') for which the photoasso- 
ciative coupling vanishes completely. This implies that 
the Jaynes-Cummings Hamiltonian produces a material 
Schrodinger cat like state of the two species, 

|*) = a|ae # ) ® \B) + b\a) ® \B'), (6) 

very analogous to the production of a similar state in 
cavity QED experiments J5j. By driving the transition 
between the active and the passive level in the B atom, 
we can subsequently monitor the interference of the two 
coherent components. The production of a mesoscopic 
superposition state can thus be verified experimentally. 
In a gauge invariant formulation of the problem, there is 
of course not a coherent state and mean atomic field in a 
single site, but a total number state of the entire atomic 
ensemble may be well approximated by a coherent state, 
and when this is split into a product of coherent states 
on each lattice site, a total number projection operator 
produces a very entangled state of the atomic compo- 
nents on each site j[7), but physical observables will be 
correctly handled by the coherent state ansatz . 

The cavity QED implementation of the Jaynes- 
Cummings Hamiltonian allows for classical excitation of 
both the two- level system and the oscillator, and it op- 
erates experimentally with the injection of a beam of ex- 
cited two-level atoms, who can give off their excitation 
energy to the oscillator degree of freedom. In the ion trap 
implementation, both the oscillator and the two-level sys- 
tem can be excited independently, and by changing the 
coupling frequency an "anti" Jaynes-Cummings Hamil- 
tonian x a<J + X* a ^ (J ^ an d a "two-phonon" Hamiltonian 
\a 2 + X*( a<! ) 2 carL be implemented to produce various 
non-classical states Q. Tomographic reconstruction of 
the motion of a trapped ion has been demonstrated in || . 
That measurement scheme made use of a classical driv- 
ing to displace the oscillator by the complex amplitude a 
and subsequently the populations of the various oscillator 
states were probed by the two-level population oscilla- 
tions (HH) as described above. The position-momentum 
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Wigner function was then obtained by an analytical for- 
mula {I from the populations recorded with different dis- 
placement a. 

It will add an experimental level of difficulty to cou- 
ple atoms coherently from a reservoir condensate into 
the lattice potentials, but with the current progress in 
transfer and loading of degenerate gasses into microfab- 
ricated structures some possibilities are certainly 

open. Without continuous loading of the lattice, there 
is another possibility to mimic the quantum optical ef- 
fects driven by external, classical fields. We illustrate 
this with a protocol for full quantum state tomography 
of the matter wave field of the A atoms: Since we can- 
not accomplish coherent displacements of the atomic field 
we suggest to use A atoms with two internal states, the 
active state |^4) and a passive state \A') which is not af- 
fected by the photoassication process. We can then use 
the \A') population as a 'local oscillator', and as long 
as the A'-population is large, coherent transfer of atoms 
between the two atomic states is well described by a co- 
herent displacement of the A atom oscillator. In fact, this 
internal state coupling is correctly described by means of 
a rotation of the Dicke collective spin of the atoms, and 
the detection of the population statistics of the |A)-state 
atoms by means of the Jaynes-Cummings dynamics can 
be used for perfect collective spin tomography: If there 
is a total of N atoms populating states |^4) and \A') in 
every lattice site, it takes AN + 1 effective spin rotations 
followed by population measurements to reconstruct the 
full collective spin density matrix, following the protocol 
applied to large hyperfine manifolds in Ref . p2| . 

In summary, we have suggested how to realize a purely 
atomic Jaynes-Cummings Hamiltonian and how to use 
it both to generate quantum optical states of matter 
and to study various many-body properties of the sys- 
tem. To isolate the different effects and mechanisms, 
we did not present studies of the joint effect of interac- 
tions, inter- well tunneling and the photoassociation pro- 
cess. It is clear that these effects in combination may 
provide wide possibilities for detection and production 
of quantum states. It has been estimated in [n| 
that one can reliably transfer atoms coherently between 



multi-mode and multi-level generalizations of the Jaynes- 
Cummings Hamiltonian. The dynamics is not restricted 
to lattices, and as an extension to the work we imagine 
that a collection of distinguishable, e.g., spatially local- 
ized B atoms immersed in a larger trapped condensates 
of A atoms can serve to induce a host of interesting ef- 
fects in analogy with the more general interaction of light 
and two-level atoms. An appealing possibility would, for 
example, be to look for a matter wave analogue of the 
Mollow fluorescence triplet, when an atomic laser beam 
interacts with an ensemble of B atoms and molecules can 
dissociate to atomic states with higher and lower energy 
than the energy of the incident beam. 
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